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OO ' The co-existence of band Jahn- Teller (BJT) effect with superconductivity (SC) is studied for 

correlated systems, with orbitally degenerate bands using a simple model. The Hubbard model for 
a doubly degenerate orbital with the on-site intraorbital Coulomb repulsion treated in the slave 
' I boson formalism and the interorbital Coulomb repulsion treated in the Hartree-Fock mean field 

f*^ . approximation, describes the correlated system. The model further incorporates the BJT interaction 

D ' and a pairing term to account for the lattice distortion and superconductivity respectively. It is 

found that structural distortion tends to suppress superconductivity and when SC sets in at low 
\^ • temperatures, the growth of the lattice distortion is arrested. The phase diagram comprising of the 

SC and structural transition temperatures Tc and Ts versus the dopant concentration S reveals that 
( , the highest obtainable Tc for an optimum doping is limited by structural transition. The dependence 

^ ' of the occupation probabilities of the different bands as well as the density of states (DOS) in the 

04 ' distorted-superconducting phase, on electron correlation has been discussed. 

P~ ! PACS. No. : 74.25.DW ; 74.20.Fg ; 76.62.-c 

■ The dependence of Tc on structural properties is rather less understood and had been a crucial question since the 
discovery of the high temperature superconductivity (HTSC) in the cuprates. On the other hand, the dominant effect 
QQ of electron correlation on the normal and SC properties of the cuprates had been realized soon after the discovery [Q. 
^\ [ The structural transitions associated with the crystallographic symmetry has been observed in many intermetallic 
"Y^ ■ compounds ^ as well as in the high Tc cuprates [|j and fuUerides Q. A common feature to all these systems is 
Ci I the existence of the Fermi level (FL) in an orbitally degenerate band. The essence of the Band Jahn Teller (BJT) 
distortion is the following. If the Fermi level of a system lies in an orbitally degenerate band then there can be a 
net gain in electronic energy if the degeneracy is removed by the lowering of the lattice symmetry due to production 
spontaneous strain which results in a structural transition. The production of strain will cost elastic energy, therefore, 
the strain will be stabilized if the gain in the electronic energy overcomes the cost in elastic energy. The system 
^ , will like to optimize the distortion so as to maximize the gain in energy, which will happen when the lower band is 
maximally occupied and the upper split band is minimally occupied. The system would however remain metallic in 
case of a small but finite distortion, as the two splited bands can overlap. On the other hand, a split band with a gap 
and a resultant insulating phase would arise because of large distortion. On deviating from the half-filled situation 
?H ' by doping holes into the system will result in lowering of the band filling, thereby reducing the gain in energy, 
hence the magnitude of strain is expected to reduce. Thus the difference in the occupation probability of the two 
bands will couple to the lattice strain which defines the Jahn- Teller Hamiltonian. However, the redistribution of the 
density of state associated with Jahn- Teller distortion is expected to affect the occurance of superconductivity in the 
system. A drastic reduction in the density of states at the FL due to BJT distortion would suppress superconducting 
transiotion. On the other hand, if the superconducting pairing is phonon mediated one can expect an enhancement 
in the transition temperature because the distortion of the lattice will soften the frequency of some phonons resulting 
in an enhancement of the coupling conatant. In this paper therefore, we present a model study of the coexistence of 
BJT distortion and superconductivity in narrow band systems. 

The electron correlation is known to play a significant role in shaping the physical properties of the normal and SC 
states of the correlated systems. It is believed that the single band Hubbard model contains all the necessary physical 
ingredients to describe electron correlation in narrow band systems. One of the recently contrived methods to solve 
the Hubbard model for any filling and any value of correlation (U) is the slave - boson technique, as formulated by 
Kotliar and Ruckenstein (KR) |5[] which introduces four auxiliary boson fields corresponding to the occupancy of a 
site, namely empty (e^), doubly occupied (d;), singly occupied {pia-/ Pi-a) with spin icr respectively. In terms of the 
slave boson field operators the single band Hubbard model takes the form 
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~4d, - pl^Pia) + J2 -^''^ ("^iTCicT - 4d, ~ pl^Pia) (1) 
(T ia 

where Zia = (1 — d|di — Pi„Pia)^^^'^{elPi<^ ~ ~ Pi-aPi-<^)~^^'^ the boson field operators are not 

independent of each other but are constrained by the requirements of completeness and local charge conservation at 
a site, hence the fourth and fifth terms are added in the Hamiltonian with the help of Lagrange multipliers (A^ and 
A^). The values of the boson field operators and Lagrange multipliers are determined by minimizing the free energy 
in the saddle point approximation. This approximation diagonalizes the U term whereas renornializes the hopping 
term as qtij with q = z ; q usually being a complicated function of the reduced correlation u {— ■^), and the 
dopant concentration (5). In this approach solutions are obtained for the paramagnetic state for all values of u and 
band fillings that reproduces the correct Brinkmann-Rice result for metal-insulator transition at a critical value of 
correlation (Uc) at half-filling. Approximately, in the weak correlation limit q = 1 — and in case of strong correlation 
and small values of S, q = _^ (for details please see ref. [^,^). 

We generalize the above formalism to a model system in which the electrons are in a two-fold degenerate Cg band 
and interact with the lattice as well as between themselves, the later via a BCS kind of pairing interaction mediated 
by some boson exchange. The Hamiltonian for such a system can be written as, 

H= ^ [{gefe -^ + 2;/' < 71^ >}nfe^„ - (-l)"Genfc^c] + ^ce^ + ^ (Ac^^^cL^^^^^ + /i.e.) 

+bosonic terms (2) 

The first term corresponds to band-energy of the correlated electrons in a two-fold orbitally degenerate Cg band ; where 
[/' is the strength of the inter-orbital Coulomb correlation treated within the mean field Hartree-Fock approximation, 
and the band renormalization factor q depends on the on-site intra-band coulomb interaction U. The second term 
within the square bracket corresponds to the symmetry breaking electron - lattice interaction that lifts the degeneracy 
of the bands by differently populating them and thereby producing the lattice strain (e). The cost in lattice energy due 
to such strain is given by the second term whereas the third term represents the usual BCS Cooper pairing between 
the quasiparticles of the correlated systems in which only intra-band pairing is considered for simplicity assuming the 
strength of pairing interactions to be of equal magnitudes . The bosonic terms in the Hamiltonian being C-numbers 
do not contribute to the dynamics of the system. In the normal undistorted phase, it is possible to show that the strain 
is exactly proportional to the population difference between the shifted sub-bands. The self-consistent equations for 
the strain (e) and the SC-order parameter (A) are given below, 

(-1) = —-tanh^^ (3) 



k.a k.a 

where Saik) — qe^ — + 2t/' < rifj > — (— l)"G'e is the band energy of the correlated normal state in the distorted 
phase and Ea{k) = \/ s'^{k) + A^ is the SC quasiparticle energy in the band a. 

In order to study the mutual influence of lattice distortion and superconductivity as a function of band filling one has 
to solve the above two coupled equations (3-4) self-consistently together with the equation for dopant concentration 
which can be obtained from the electron number conservation (2 — (5) = E] < "^fecra >• In order to solve these three 

kaa 

equations the sum over k as usual is converted into an integral over energy variable e, with an approximate peaked 
density of states No{e) at the centre of the unperturbed bands as given below. 
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N,ie) = NiO)Jl-\- \ In I — I (5) 



where A^(0) is the unperturbed density of states of the free electron system and 2B (which is chosen as 1 eV) is the 
band width. 

Figure 1. shows the effect of electron correlation on the temperature dependence of the strain, i.e, e(T). In order to 
see this the curves for e(T) corresponding to the value of the coupling constants G = 1.2 and band occupation n = 0.9 
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are plotted for different values of u and u' . There is a suppression of strain around 40 K due to the appearance of 
superconductivity when u and u' are zero. As the intra-orbital intra-site Coulomb correlation is switched on with 
u — 0.2, while keeping u' = 0, the strain as well as Tg has increased, and the SC transition temperature is lowered to 
around 22 K. Furthermore, when the intcr-orbital intra-site Coulomb interaction strength is increased from zero to 
u' = 0.001, there is a further suppression of Tc and an enhancement of as well as strain as can be seen from the 
dash - dotted curve. On further increasing u' to 0.005 superconductivity almost vanishes (dotted curve) while Tg and 
strain increases by a larger magnitude. Therefore, the present results predict that the suppression of strain will be 
lessened with increasing correlation in the superconducting state. 

An unusual behavior in the temperature dependence of the SC order parameter A(T) and strain e(T) is observed 
when the value of the SC transition temperature in the absence of strain [T^] is comparable to the structural transition 
temperature T^, which is depicted in Fig 2. The figure depicts the temperature dependence of e and A for the same 
value of the coupling constant G=1.2eV as earlier but for a band filling n = 0.7. In this case the strain appears at 
the temperature of Tg « lAQK, increases with decreasing temperature but drops precipitously to zero around 92 K, 
as soon as superconductivity appears. The SC order parameter rapidly increases within an extremely narrow range 
of temperature but it vanishes again with the recurrence of strain which attains a large value when the temperature 
is reduced slightly. Such alterations of SC and strain persist until the temperature is slightly less than 60 K below 
which the superconducting state stabilizes and the strain vanishes. This corresponds to a metastable situation that 
persists within a narrow range of temperature. While the SC state stabilizes at lower temperatures, the strain exists 
in the higher temperature range, with the two never co-existing in their respective regions of stability. 

In figure 3, the occupation probability/per spin for different orbitals are plotted as a function of temperature. The 
occupation probability of the two orbitals are the same above the structural transition temperature Tg but below , 
the occupation probabilities being different creates a population difference between the orbitals and hence a net strain 
builds up. With lowering of temperature, the system undergoes SC transition and hence the electron density in each 
orbital is suppressed below Tc- Needless to say, this transition resembles to the orthorhombic to tetragonal transition 
in high- Tc systems. However, as a result of reduction in the orbital population density below Tc, the growth of the 
spontaneous strain gets arrested. 

In figure 4, the detailed phase diagram comprising of Tc & Ts as a function of hole concentration is shown for 
different values of intra and inter orbital correlation strengths. The phase diagram (solid lines represent Tc curves) 
demonstrates that the Tc is largely suppressed in presence of structural distortion and the highest Tc is obtainable for 
an optimum doping ((5c) where Tg vanishes. It is to be noted that the suppression in Tc is maximum at half-filling as 
the strain is the largest there. These results are in close agreement with experimental observations [||J^. 

The structural transition in this model results from the competition between lowering of the electronic energy due 
to lifting of the orbital degeneracy (band Jahn- Teller effect) and the increase in elastic energy associated with the 
apperance of strain. The gain in energy in this process stabilizes the structural transition and this gain is largest 
when the Fermi level (FL) in the undistorted phase lies at the center of a peak in the DOS, because in that case a 
large number of electrons in higher energy states are transferred to lower energy states. This is realized in the present 
model at half filling of the degenerate band resulting in a splitting due to JT effect by an amount 2 Ge which leaves the 
lower band almost fully occupied whereas the upper band is nearly empty. Thereby, it creates maximum population 
difference between the two shifted sub-bands resulting in the largest strain. Consequently, the DOS at the FL gets 
reduced drastically so that the system is no longer favorable for superconductivity. However, as one moves away from 
half-filling of the bands (by doping holes) the FL moves away from the peak in the DOS resulting in a lesser gain in 
electronic energy which will tend to reduce the distortion as well as the structural transition temperature. As a result 
the DOS at the FL will increase thereby favoring superconductivity. Thus the SC-Tc increases with increasing hole 
concentration. 

On the other hand, the appearance of SC reduces the occupation probability of the bands exponentially due to the 
presence of the gap 2A in the energy spectrum resulting in the reduction of the strain (e) below Tc. This explains the 
arresting of the growth of the strain at Tc and its suppression below Tc. Such a picture becomes self-evident from Fig. 
3. Our results presented in figures have close resemblence to experimental data of high temperature superconductors 

Finally, within the present formulation the electron correlation has the following effects. The band narrowing due 
to renormalizatin factor q enhances the density of states which favors the Jahn - Teller splitting, while the pairing 
amplitude gets suppressed, due to the reduction of the strength of the interaction by a factor if'. Inclusion of inter 
orbital coulomb correlation makes the J-T splitting asymmetric with respect to the centre of mass of the original 
band and hence enhances the lattice strain but suppresses SC further. Therefore, the present model study helps in 
understanding the interesting phenomenon of the interplay of structural transition with superconductivity in corre- 
lated systems like the cuprates. However, it remains to be seen how the antiferromagnetic spin fluctuations leading to 
d-wave pairing superconductivity is influenced by structural distortion. This is the subject matter for a future study. 
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Figure Captions 

Fig. 1 The thermal variation of the lattice strain (e) for n = 0.9 and G = 1.2 in presence of on-site coulomb correla- 
tion (intraband (w) as well as interband {u' is denoted as ul in figures)). 

Fig. 2 The temperature variation of the superconducting gap and the stain for n = 0.7, for G = 1.2. A mctastablc 
thermal regime from around 60 K - 90 K is worth noticing (the solid line indicates the Agc{T) and the dotted one 
e(T). 

Fig. 3 The temperature variation of the occupation probabilities of the splitted bands for n = 0.9, G = 1.2 in 

presence of the onsite intraband (u) as well as (ul = u') coulomb correlations. 

Fig. 4 The complete phase diagram in terms of the structural and superconducting transition temperatures (Tg, Tc) 
as a function of hole concentration (6) for G = 1.16. The curves with connected lines correspond to Tc curves). 
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